function [pop, y, best] = initializeCC(NP, N, upperBound, lowerBound, fun)
pop = lhsdesign(NP, N) .* (upperBound - lowerBound) + lowerBound .* ones(NP, N);
y = fun(pop);
[~, ind] = min(y);
best = pop(ind, :);
end

